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We present the geodesies on homogeneous and isotropic negatively curved spaces in a simple form suitable for application 
to cosmological problems. We discuss how the patterns in the microwave sky of anisotropic homogeneous universes can be 
predicted qualitatively by looking at the invariances that generate their three-geometries. 



ON . 

o\ ■ 

Oh- 

<: 

m : 



o 
o 

On 



Of 
i 1 

O 

CO 



X 



98.70 Vc, 98.80.Cq, 98.80.Hw 



Metric theories of gravity, of which Einstein's general 
theory of relativity is the most elegant known example, 
couple the material content of the Universe to its space- 
time geometry. As a result the paths of light rays reflect 
the curvature of space and provide a means to deter- 
mine the density and overall curvature of the Universe. 
Locally, the effect of the Sun's gravitational field upon 
grazing light rays and radio waves provides a powerful 
test of Einstein's predictions. Globally, we have begun 
to observe the effects of gravitational lensing, and we 
know that varying the density of the Universe will al- 
ter the observed features of microwave background pat- 
terns traced by photons on the sky. Although most 
versions of the inflationary universe scenario lead us to 
expect that the universe will be expanding very close 
to the critical divide, with the density parameter sat- 
isfying | fio — 1| ~ O(10 -5 ), there are varieties of infla- 
tion which predict that the universe is significantly open 
(f2o < O(0. 1)) jjj. Moreover, the observational evidence 
stubbornly refuses to provide a clear endorsement of the 
Oo ~ 1 predictions. Primordial nucleosynthesis limits 
the baryon density to fall well short of this value, and 
so non-baryonic forms of dark matter must be found in 
support of a closure density. For these reasons, the de- 
termination of all the observational differences between 
flat, closed, and open universes is an important goal for 
cosmologists. The most sensitive discriminators promise 
to come from the study of null geodesies. 

There has been much recent interest in studying the 
observational signatures of open universes that possess 
compact topologies which are produced by the periodic 
identification of space coordinates, and of non-compact 
topologies that possess identifications in some space di- 
rections which permits them to have integrable 
geodesic motions. Open universes are interesting can- 
didates for possessing non-trivial topologies because the 
curvature of space provides a natural length scale to re- 
late to the scale of topological identification. There have 
been several investigations of the effects of topological 
identifications on the power spectrum of the microwave 



sky and on simulated COBE sky maps. These studies 
enable us to limit the scale of topological identifications 
more powerfully than by searching for multiple images of 
prominent luminous sources In these studies, the neg- 
ative curvature of open universes plays an important role. 
Geodesic flows on compact negatively-curved spaces are 
chaotic, as first noted by Hadamard ||, and the geodesic 
flow on a compact negatively-curved space has become 
a key paradigm for the identification of chaotic classical 
motions as well as in the study of quantum chaos. Some 
possible implications for the microwave background have 
been discussed in ref. However, the formal character- 
isation of these mixing flows ]?]] is of little practical utility 
for the study of geodesies in the universe. The needs of 
cosmologists are more specific and in this paper we aim to 
provide the geodesies for open universes in simple usable 
form. These do not seem to exist in the literature and can 
be employed to study the behaviour of the cosmic back- 
ground radiation in open universes with any topology. 
With such applications in mind, we provide a brief anal- 
ysis of the geodesies on an expanding negatively-curved 
space. They are presented in an immediately accessible 
form for use in cosmological investigations. Finally, we 
provide some new qualitative discussion of the microwave 
sky patterns in anisotropic universes. 



II. GEODESICS ON A STATIC, NEGATIVELY 
CURVED SPACE 



The most familiar description of the metric in an open 
universe is in hyperbolic coordinates 



ds 2 



-dt 2 + dr 2 + sinh 2 r (d6 2 + sin 2 



(1) 



However, the geodesies can be found more simply in a 
coordinate system (x, y, x) related to (r, 9, 4>) by 



e~ z = cosh r — sinh r cos ( 
e~ z x = sin cos <fi sinh r 
e~ z y = sin 6 sin 4> sinh r . 



(2) 
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In the (x, y, z) coordinate system the metric is 



W 2 = z = Zi ± 777, 



(12) 



-dt 2 + dz 2 + e- 2z (dx 2 + dy 2 ) 



(3) 



The geodesic equations can be found in the usual way 
from 

but it is more efficient to introduce the Lagrangian 

L= l -[-t' 2 + z' 2 + e- 2z (x' 2 +y 12 )] (5) 

where ' = d/dX. The equations of motion can be found 
from 



, dL 



(6) 



where IT g = dL/dq' is the momentum conjugate to coor- 
dinate q. 

Since the Lagrangian is independent of (t,x,y), the 
corresponding conjugate momenta are conserved, giving 



n t = -t' = -Ei 
n, = e- 2z y' = U yi 



(7) 



where Ei,Tl x i,Ilyi are constants of the motion. The last 
of these equations is 



// 1 — 2zi 12 1 I2\ n 

z + e [x + y ) = 



(8) 



This second-order equation can be reduced to first or- 
der by exploiting the invariancc of the length element. 
Choosing the affine parameter to be proper time 77, 
ds 2 /di] 2 = a (a — for photons and a = — 1 for massive 
particles), we have ([|) 



a 



-E 2 + z' 2 + e- 2z (x' 2 +y' 2 ) 



(9) 



Substituting the solutions for x 1 and y' and defining 
W 2 = n^j -(- n^; , we can solve this for z' to find 



2z w 21 1/2 



z' = ± [(a + £?) - e 2z Wi 2 ] 



(10) 



The system has been reduced to first-order equations. 
Notice that z" < always. It follows that if z' > 0, then 
z will reach a maximum and then the geodesic reverses 
direction. 

Integrating the z equation, we have 



dz 



zi [{a + E 2 )-e^W 2 ] 1/2 



Hv-m)- (11) 



This completes the solution. Let the initial time be rji 
0. Firstly, if W 2 — 0, the trajectories are simple lines, 



with 7 = (a + E 2 ) 1 / 2 . If Wi ^ we can integrate (p 
and then (0) to obtain 



W 2 e 2z = 7 2 



1 



cosh (7?7=F/3i) 

wi 



2 7 [tanh(7?7 =p ± tanh(/?i)] 



y = W + [tanh( 7 77 T A) ± tanh(A)] . (13) 

The last constant of integration that appears here is 

tanhA = \j 1 2 - e 2zl W 2 h ■ (14) 

Notice that, unless a particle is shot directly along the 
axis, it will never reach z — +00. The constant, which is 
proportional to tanh(A), could be absorbed into X{, but 
then Xi and would not correspond to the value x(r)i = 
0). Rewriting these in another way, that is useful when 
tracing geodesies, gives 



z = e Zi 7Cosh(7?7) =p smh.{ , yri)'J~^( 

1 TT 2zi 2 

x = Xi + ll^ie '7 



tanh(7?7) 



= yi + Il yi e 2z ^ 2 



where 
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7 =p tanh(777)-y/7 2 — e 2zi W 2 

tanh(7?7) 
7 =p tanh(77y)A/7 2 - e 2Zi W 2 



exp(2z;) 



cosh 2 /3iWj 2 
We can write the path parametrically as 



s a * = ^-(«-*i) 2 -(l/-w) a 



(15) 



(16) 



(17) 



where x- t = X; ± (n a .i/W i 2 )7tanh(/3i) and is defined 
analogously. 

These trajectories have some very odd properties. As 
already mentioned, the trajectories never reach z = +00. 
Any photon travelling along increasing z eventually hits 
a maximum and then wraps back. As it moves toward 
z = — 00, the velocities along x and y fall (eqns (0)) so 
that they never reach infinite values of x or y. 



A. Geodesic flows in an expanding universe 

When the expansion of space is included, the full met- 
ric becomes that of the open Friedmann universe, 

ds 2 = -dt 2 + a 2 (t) [dz 2 + e- 2z (dx 2 + dy 2 )] . (18) 
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The null geodesic equations are 



III. TRACING GEODESICS 



t + HF =0 



+ e~ 2z (x 2 + y 2 ) + 2Hz = 



x - 2zx + 2Hx = 



y - 2zy + 2Hy = 



(19) 
(20) 
(21) 
(22) 



where an overdot denotes d/dX and Ti — Ht = dlna/dX. 
Notice from (Il9|) that 



t = 



1 



(23) 



so dX = adt. 

All of the results of the previous section can quickly be 
adapted to the case with expansion if a time coordinate 
is chosen astutely. Now let / = d/d-q. Then the geodesic 
equations become 

rf \z" + e- 2z (x' 2 + y' 2 )} + z' [ij + 2Hfj\ = (24) 



fj 2 [x" - 2y'x'\ + x' [fi + 2W]\ = 

ri 2 [y" -2z'y']+y'[f ] + 2W 1 ] =0 . 
If we choose the coordinate r\ such that 

r\ + 2Ui] = 0, 
then it follows that the geodesic equations become 
z" + er 2z (x' 2 + y' 2 ) = 



(25) 



(26) 



(27) 



V 

t" - 2z'x' = 
,,// 



(28) 



y 



2z'y' = 



which is precisely the same as those on a static, negatively 
curved hypersurface. The solutions are then the same 
as equations (|l^) (or equivalently (|l5|)) with 7 = 1 for 
photons and 



V 



dX 



dt 
W) 



(29) 



is the usual conformal time. 

For timclikc geodesies, the motion can again be pro- 
jected onto a static hypersurface with a new affine pa- 
rameter. Again, we recover the geodesic flows of the pre- 
vious section with a — — 1 except the time parameter for 
massive particles is not conformal time but rather m 



V = 



dt 



i(v + a 2 ) 1 / 2 



(30) 



and v a = av{t)/\/\ — v 2 and v 2 (t) — giji l xK 



The microwave background provides a sensitive probe 
of the curvature of the universe. If we locate the Earth (or 
our near-Earth satellite) at the origin of the coordinate 
system, then only photons from the surface of last scat- 
ter which travel along radial geodesies will be observed. 
The geodesic equations are then greatly simplified with 
respect to the direction of observation on the sky. 

The photons seen in the sky today can be traced back- 
wards to locate their point of origin. We have six un- 
knowns: 



Zi < = > Pi 
Xi 



(31) 



The magnitude of z[ is fixed by eqn (|0|) once the other 
initial coordinates are specified. The choice z[ corre- 
sponds then to a choice of the sign in eqn ([l(]). There 
are two sets of boundary conditions: (i) The photon po- 
sition vector is at the origin, and (ii) The velocity vector 
is opposite to the unit vector pointing in the direction of 
observation. In other words, 



(i) x = x a = 
(ii) v — — h(Q, 4>) 



(32) 



Boundary condition (i) places the Earth at the origin of 
the coordinate system. Using the geodesic solutions, we 
have 



cosh (r\ =F 0i) 



Xi = - 



Ih 



[tanh(ry Q =p #) ± tanh(/?i)] 



IV 

^ [tanh(ry =F Pi) ± tanh(A)] 



(33) 



When evaluated at the origin, the geodesic equations re- 
late the components of the velocity vector today to their 
initial values: 



12 = 1 _ w 2 

o ^ vy 1 



Vo = U yi 



(34) 



This velocity vector is normalized to 1 as it must be for 
photons. Using 



sin 6 cos 4>\ + sin 6 sin 0j + cos 0k , 



(35) 



we can rotate this into the (x, y, z) coordinate system at 
the origin to find 



z <=> k 
x <=> j 
y<=>l . 



(36) 



3 



Condition (ii) then gives 

z'= — cos 9 



x' Q = U xs = — sin 9 cos ( 
y' = Uyi = -sin6»sin<; 



(37) 



> It follows that W 2 = sin 2 9. Also, taking the derivative 
of the geodesic equation, 



z' Q = — tanh(?7 qp A) = — cos#, 

from which it follows that 

cos 9 — tanh(?/ ) 



=Ftanh(/3i) 



1 — tanh(?y ) cos 9 



(38) 



(39) 



Pi = ±(?7 — arctanh(cos6>)). Putting ( |37| ) into ( p3| ) we 
find 



e~ Zi = cosh(?7 ) — sinh(r7 ) cos( 
e~ Zi Xi — sin0cos0sinh(ry o ) 
e~ Zi yi — sin0sin0sinh(?7 o ) . 



(40) 



These are all radial geodesies in (r,0,<fr); that is, r = 
rj — r). The integrated Sachs- Wolfe effect considers only 
radial geodesies; but processes such as gravitational lens- 
ing, that can deflect a photon into the line of sight, would 
draw from the more general pool of non-radially directed 
photons. 



IV. LOCAL AND GLOBAL ANISOTROPY 

The full geodesies on a universe of negative curvature 
have been obtained in an explicit form most accessible to 
cosmologists. One arena of renewed interest where the 
full geodesies may be needed is the case of a small uni- 
verse. Negatively curved spacetimes can be made small 
and finite through topological identifications. A small 
universe could be witnessed with periodic effects or by 
features in the power spectrum of the microwave back- 
ground. Topology induces global anisotropy even when 
the underlying space is locally isotropic. Local anisotropy 
is also possible in the absence of topological identifica- 
tions. When homogeneous anisotropies are present, ei- 
ther in the form of shear or rotation in the expansion 
of the universe, or possibly also in the three-curvature 
of space, there are only a finite number of homogeneous 
spaces which can provide an exact description of the ge- 
ometry of space. These anisotropic spaces were first clas- 
sified by Bianchi H. They were introduced into cos- 
mology by Taub 110], and presented in the most effi- 
cient manner by Ellis and MacCallum pj. Since the 
microwave sky is currently the most significant histori- 
cal record of the primordial radiation, it is instructive to 
show how the anisotropic sky patterns created by these 
different anisotropic universes can be predicted just from 



a knowledge of the group invariances that generate the 
homogeneous geometries and their geodesic flows. 

In order to determine the detailed sky patterns per- 
mitted by the Bianchi geometries it is necessary to solve 
for the evolution of the geodesies on the anisotropic cos- 
mological models either exactly or approximately (in 
the case of small anisotropy); see for example refs. 
Jl2|,^3 14,16|fL7||. Again, the most unusual features arise 



in open universes. The basic quadrupole pattern arises in 
the simplest (flat) Bianchi type I universe with zero cur- 
vature. The addition of negative curvature focuses this 
quadrupole into a small hotspot on the sky (there is a pre- 
ferred direction because there is a direction of lowest and 
highest expansion rate) in type V universes, which still 
possess isotropic 3-curvature. If anisotropic curvature is 
added then we reach the most general class of anisotropic 
homogeneous spaces and a spiralling of the geodesies is 
added to the quadrupole or focused quadrupole in the 
flat or open type VII universes. These particular models 
have been studied in the past by linearizing the geodesic 
equations about the isotropic solutions in which the tem- 
perature anisotropy of the microwave background is zero 



3 14 



,|bj,|17j,|18|,yj|. They describe the most gen- 
eral anisotropic distortions of flat and open Friedmann 
universes. However, it is also possible to predict the ge- 
ometric sky patterns expected in the different Bianchi 
type universes by simply noting the nature of the groups 
of motions which define each homogeneous space. 

The Bianchi classification of spatially homogeneous 
anisotropic universes is based on the geometric classi- 
fication of 3-parameter Lie groups. The action of these 
groups on the spacelike hypersurfaces of constant time 
in these universes can be prescribed by three transfor- 
mations of cartesian coordinates [x, y, z). Each model 
possesses two simple translations in the x — y plane, with 
generators d/dx and d/dy, together with a more compli- 
cated motion out of this plane which is different for each 
group type. If it is considered as a flow from the z — 
plane to some other plane, z = a = constant, then the na- 
ture of this flow tells us qualitatively what the microwave 
background anisotropy pattern will look like. In the sim- 
plest flat universe of Bianchi type I the z— flow is uniform 
and just maps (x, y) — > (x, y). This corresponds to a pure 
quadrupole geodesic temperature anisotropy pattern. In 
the open Bianchi type V universe the z— flow is a pure di- 
lation and maps (x,y) — > e a (x 1 y), with a constant. This 
dilation describes the hotspot created by the focusing 
of the quadrupole pattern in open anisotropic universes. 
The most general non-compact homogeneous universes 
containing the Friedmann models are of Bianchi type 
VII/j, which contain the flat Friedmann universes when 
h — and the open Friedmann universes when h 0. In 
type VII/j the z— flow is a rotation plus a dilation: circles 
of radius r = (x 2 + y 2 ) 1/>2 = 1 are mapped into circles of 



radius r 



,avh 



and rotated by a constant angle a. Ob- 
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servationally, this corresponds to the geodesies producing 
a focusing of the basic quadrupole (as in type V) with a 
superimposed spiral twist. In Bianchi type VIIo there is 
simply a spiral added to the underlying quadrupole with 
no focusing because the 3-geometry is flat. The z— flow 
for the closed universes of Bianchi type IX, which con- 
tain the closed Friedmann models as special cases, is more 
complicated. The z— plane action corresponds to the fol- 
lowing 50(3) invariant motion in polar coordinates based 
on the x-axis, in the region r < 2n, [ pp[ : 

{r,6, 4>} = {cos- 1 ^ - /3 2 ) 1 / 2 cos(a/2)], 

sin -1 R[l - (1 - 1 ) cos 2 {a/2)}- 1 /a j i a + O } 

Sky patterns in other, less familiar, Bianchi types can be 
generated in a similar fashion, if required. Thus, type 
Ylh is generated by a z— flow that combines a shear with 
a dilation: the hyperbolae x 2 — y 2 = A 2 are mapped into 
hyperbolae which are rotated by a hyperbolic angle (3 
into x 2 - y 2 = A 2 e 2f3 ^ r \ 
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